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Using full-potential linear augmented plane wave method (FP-LAPW) and the density functional 
theory, we have carried out a systematic investigation of the electronic, magnetic, and cohesive 
properties of the chalcogenide CrTe in three competing structures: rock-salt (RS), zinc blende (ZB) 
and the NiAs-type (NA) hexagonal. Although the ground state is of NA structure, RS and ZB are 
interesting in that these fcc-based structures, which can possibly be grown on many semiconductor 
substrates, exhibit half-metallic phases above some critical values of the lattice parameter. We find 
that the NA structure is not half-metallic at its equilibrium volume, while both ZB and RS structures 
are. The RS structure is more stable than the ZB, with an energy that is lower by 0.25 eV/atom. 
While confirming previous results on the half-metallic phase in ZB structure, we provide hitherto 
unreported results on the half-metallic RS phase, with a gap in the minority channel and a magnetic 
moment of 4.0 fiB per formula unit. A comparison of total energies for the ferromagnetic (FM), non- 
magnetic (NM), and antiferromagnetic (AFM) configurations shows the lowest energy configuration 
to be FM for CrTe in all the three structures. The FP-LAPW calculations are supplemented by 
linear muffin-tin orbital (LMTO) calculations using both local density approximation (LDA) and 
LDA+U method. The exchange interactions and the Curie temperatures calculated via the linear 
response method in ZB and RS CrTe are compared over a wide range of the lattice parameter. The 
calculated Curie temperatures for the RS phase are consistently higher than those for the ZB phase. 

PACS numbers: 71.20.Nr, 71.20.Lp, 75.30.Et, 75.10.Hk 



I. INTRODUCTION 



Half-metallic (HM) ferromagnetism^ is an ideal prop- 
erty of materials from the viewpoint of their applica- 
tion in spintronic devices.— — Thus a considerable re- 
search has been devoted in recent years to the design- 
ing and measurement of the physical properties of such 
materials. One group of such materials is known as 
the dilute magnetic semiconductors (DMS)^— where the 
magnetic properties are mediated by the percolation ef- 
fects connecting a very small concentration of transition 
metal (TM) atoms such as Mn or Cr in a semiconduc- 
tor host. Another group, that has been studied exten- 
sively, consists of TM-semiconductor binary/ternary al- 
loys, where the TM concentration can be as high as the 
other components^— In particular, some ZB structure 
transition metal chalcogenides and pnictides have been 
predicted to be HM ferromagnets on the basis of theoret- 
ical calculations*^— Here the magnetic properties of the 
transition metal matrix is modified via the hybridization 
between the transition metal 3-d and the s and p orbitals 
of the pinctogen or the chalcogen. Half-metallicity arises 
as a result of interplay between the exchange splitting 
and the hybridization. Some Cr-based chalcogenides and 
pnictideS ) 17 i 24 ~ — have been predicted to be not only HM 
but also potential candidates for possessing high Curie 



temperatures T c . Akinaga et al~ have reported grow- 
ing ZB thin films of CrAs on GaAs (001) substrates via 
molecular beam epitaxy. These showed ferromagnetic be- 
havior at temperatures in excess of 400 K and magnetic 
moments of 3/x# per CrAs unit. Thin films of CrSb grown 
similarly on GaAs, (Al,Ga)Sb, and GaSb substrates have 
been found to be of ZB structure and ferromagnetic with 
T c higher than 400 K. 35 Deng et aL— were successful in 
increasing the thickness of ZB-CrSb films to ~ 3 nm using 
(fn,Ga)As buffer layers, and Li et al~ were able to grow 
~ 4 nm thick ZB-CrSb films on NaCl (100) substrates. 
Two of the present authors have recently presented an 
extensive theoretical study of the exchange interactions 
and Curie temperatures in a series of Cr-based pnictides 
and chalcogenides in ZB structure as a function of the 
lattice parameter, 3 - 

One important issue is that although the ZB phase of 
these compounds shows half-metallicity and often (the 
promise of) high Curie temperature, the ground state 
usually has a different structure. For many of these com- 
pounds the ground state is NiAs-type hexagonal (NA). 
The semiconductor substrates on which thin films of 
these compounds are grown have cubic structures, usu- 
ally RS or ZB. The other point is that often the equi- 
librium volume does not show half-metallicity, which ap- 
pears only at an expanded volume. From this stand- 
point it would be useful to explore the possibility of the 



half- metallic phase in structures other than ZB, and in 
particular the RS structure. In this work we have ex- 
amined this problem for CrTe. Zhao and Zunger— ar- 
gued that ZB CrTe should be epitaxially unstable against 
the NA structure. However, recently Sreenivasan et 
alM' and Bi et al.- have reported growing ^5 nm 
thick thin films of ZB CrTe by low temperature molecu- 
lar beam epitaxy. Their measurements of temperature- 
dependent remanent-magnetization indicate a Curie tem- 
perature of ~100 K in these films. One attractive feature 
of bulk CrTe, at least from a theoretical view point, is 
that the equilibrium volume in the ZB structure is al- 
ready half-metallic, whereas most other Cr-based chalco- 
genides and pnictides show half-metallicity only at ex- 
panded volumes. 17 Besides, the half-metallic gap is larger 
in the chalcogenides like CrTe or CrSe than in the pnic- 
tides such as CrAs or CrSb. This is due to larger Cr- 
moment in the half-metallic phase for the chalcogenides 
(^4/xs) than for the pnictides (~3/is), which results in 
larger exchange splitting. Our previous study-2 7 . has re- 
vealed that chalcogenides such as CrS and CrSe have a 
high likelihood of exhibiting antiferromagnetic or com- 
plex magnetic behavior for a range of the lattice parame- 
ter, whereas such tendencies are less pronounced in CrTe. 
Hence, the focus of this work is on CrTe. We study the 
appearance of half-metallicity in three competing struc- 
tures for CrTe: NA, RS and ZB. The ground state struc- 
ture (NA) shows half-metallicity not at the equilibrium 
volume but at some expanded volume, while both RS and 
ZB structures show half-metallicity at their respective 
equilibrium (lowest energy) volumes. The RS structure is 
found to be more stable than the ZB, with an energy that 
is lower by 0.25 eV/atom. Some band calculations for 
CrTe in the NA structure have already appeared 4 ^ in the 
literature, and the same goes for ZB CrTej 13 ' 17 Ahmadian 
et al.— have recently discussed half-metallicity at a ZB 
CrTe(OOl) surface and its interface with ZnTe(OOl). So 
far CrTe in the RS structure has not been studied from 
the viewpoint of its electronic and magnetic properties. 
Although both RS and ZB structures are fec-based, their 
magnetic properties should differ somewhat for the same 
lattice parameter due to different levels of hybridization 
between the Cr 3-d and Te sp-orbitals. This should influ- 
ence the magnetic moments, half-metallic gaps and the 
exchange interactions. We find that the exchange inter- 
actions are predominantly ferromagnetic and strong in 
the RS structure, while in the ZB structure Cr-Cr inter- 
actions can have non-negligible antiferromagnetic com- 
ponents for some lattice parameters. This results in a 
lower Curie temperature T c for the ZB structure over a 
wide range of the lattice parameter. The promise of lower 
energy and higher T c should prompt the experimentalists 
to try to fabricate CrTe thin films in the RS structure. 



II. ELECTRONIC STRUCTURE 

The two fec-based structures ZB (F43m, No. 216) 
and RS (Fm3m, No. 225) are well-known. The NA 
structure(P63/mmc, No. 194) is based on a distorted hep 
array of the arsenide anions. By contrast with the 
wurtzite structure, which shares the same formula type 
(AB) with the NA structure, the cations occupy all the 
octahedral sites rather than half the tetrahedral holes. 
There is one octahedral hole for each hep lattice site, 
and so the AB stoichiometry is preserved. The local co- 
ordinations of the anions and cations are different in the 
NA and wurtzite structures. Bulk samples of CrTe are 
known to crystallize in this structure 4 ^ with cell param- 
eters a = 3.997A, c = 6.222A at room temperature. 
The Cr-Te bond length is 2. 78 A and the c/a ratio is 
1.557, different from the ideal hep value of 1.633. Thus 
the six Te atoms form a trigonally distorted octahedron 
around Cr. The Cr-Cr nearest neighbor distance in the 
ab plane is a = 3.997A, while the out of plane distance is 
c/2 = 3.11lA. Readers may consult Fig. 1 of Ref.[Il[ to 
view the structure. 

For the electronic structure study we have used the 
WIEN2k code , 44 ' 45 which employs the full-potential lin- 
ear augmented plane waves (FP-LAPW) plus local or- 
bitals method. The generalized gradient approximation 
(GGA) proposed by Perdew et al. was used for ex- 
change and correlation potentials^ According to a study 
by Continenza et al.^ GGA is essential for obtaining 
accurate equilibrium volumes and magnetic moments for 
Mn-based pnictides of ZB and NA structures. We con- 
sider full relativistic effects for the core states and use 
the scalar approximation for the valence states. Spin- 
orbit coupling is neglected, as it should have negligible 
effect on the results for CrTe. Convergence with re- 
spect to basis set and k-point sampling was carefully 
verified. About 3000 K points were used for the Bril- 
louin zone (BZ) integrations, using the Monkhorst-Pack 
scheme^ and 14 x 14 x 14 divisions of the BZ. We used 
R m t * K max = 8.0 and angular momentum expansion up 
to Imax — 10 in the muffin tins, and G ma x = 14 for the 
charge density. All core and valence states were treated 
self-consistently. For the calculations involving the dif- 
ferent structures the muffin-tin (MT) radii were taken 
to be equal to 2.2 and 2.5 bohrs for Cr and Te atoms, 
respectively. The same MT radii were used in the calcula- 
tions of different structures at their equilibrium volumes 
in order to compare accurately the total energies. Self- 
consistency in charge density was achieved to an accuracy 
higher than 10~ 4 . 

We carried out a systematic structural optimization 
of NA, ZB, and RS CrTe in FM, NM and AFM config- 
urations. For the NA structure the c/a ratio was op- 
timized and the result c/a = 1.542 for the FM phase 
is in good agreement with previously obtained result^ 
This optimized c/a ratio was used for the final NA struc- 
ture volume optimization. Bulk modulus was calculated 
by fitting the volume vs. energy results to the empiri- 



cal Murnaghan equation of stated 9 We find that for all 
three structures the FM phase has lower energy than the 
NM and AFM phases. To assess the relative stability of 
CrTe in different structures we consider the total energy 
per formula unit in the ZB and RS structures relative 
to the NA structure. The calculated optimized equilib- 
rium lattice constant a, bulk modulus B, the total en- 
ergy difference Ei between the NM and the FM states: 
Ei = E NM - E FM , and E 2 between the AFM and the 
FM states: E 2 = Eafm — E FM , and the metastabil- 
ity energy E m , which is the total energy per formula 
unit relative to the NA structure, are given in TableHTl 
Note that for RS and ZB structures, in calculating E 2 
we have considered the energies of three AFM configura- 
tions: [001], [110] and [111] and considered the lowest of 
these energies. For the hexagonal NA structure, the AFM 
state was created by simply ascribing equal but opposite 
magnetic moments to the two Cr atoms in the unit cell. 
Cohesive properties are dependent on the treatment of 
the exchange-correlation effects in a non-trivial manner, 
as indicated by the LDA values of the lattice parameter 
and bulk modulus shown in parentheses in TabldlTlfor the 
RS and ZB equilibrium structures, where all other values 
are from GGA. Increased correlation effects in GGA in- 
creases the equilibrium lattice parameter and lowers the 
bulk modulus. 

It is apparent from TablcHT] that the FM state is lower 
in energy than the NM and AFM states for all the three 
structures and the ground state structure is NA. Our cal- 
culated CrTe results of ZB and NA by GGA are in good 
agreement with previous results i 13 ' 42 The calculated total 
energy as a function of the atomic volume for FM phase 
ZB, RS, and NA structures of CrTe is shown in FigfT] 
The ZB and RS phases are higher in energies than NA 
by 0.328 eV and 0.076 eV/atom, respectively. However, 
the RS structure is only marginally unstable against the 
NA structure, and more stable than the ZB structure, 
with an energy that is about 0.25 eV/atom lower. Thus 
we have carried out a detailed study of the RS phase. For 
lack of both experimental and theoretical data regarding 
RS CrTe, the present results should be considered as pre- 
dictions and are intended to serve as reference for future 
experimental work. 

As the density of states (DOS) for the ZB structure 
has been described in several publications ) 13 ' 17 ! 32 ' 33 we 
present only the DOS's for the RS and NA structures in 
this work. For the RS structure the majority-spin chan- 
nel, as in case of the ZB structure, is metallic, while the 
minority-spin states have an energy gap of width 0.76 
eV around the Fermi level. The minority-spin gap for 
the ZB structure is 2.9 eV. All of these values refer to 
the respective equilibrium volumes. Thus both ZB and 
RS phases are HM systems. The HM gap, which is de- 
fined in Ref. 13[ is 0.83 eV for the equilibrium ZB volume, 
while for the RS equilibrium volume it is very small, be- 
ing about 0.03 eV. However, it increases with increas- 
ing volume. The spin-polarization at the Fermi level is 
quite small for RS CrTe, but the half-metallicity is robust 
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FIG. 1: (Color online) Energy as a function of volume per 
atom in the three structural phases of CrTe. 



against volume expansion. The general structure of DOS 
is similar around the Fermi level for both structures. A 
closer inspection, however, reveals some differences. In 
the spin-up channel there is a gap in the energy interval 
between 1.62 and 2.82 eV just below the Fermi level for 
the ZB structure. On the other hand, for the spin-down 
DOS the gap at the Fermi level in ZB CrTe is larger 
than that in RS CrTe. Also the peak structures differ 
by an energy shift with respect to the Fermi level. All 
these differences should be accountable on the basis of 
different levels of hybridization between the Cr-3<i and 
Te-5p orbitals and different exchange splittings for the 
two structures. The growth of metastable ZB-CrAs, 24 
ZB _ CrSb 35 i 36 and ZB-CrTe 3 ^ indicates the suitability 
of RS-CrTe for growth on carefully chosen substrates. 
The equilibrium lattice constant of RS CrTe (5.727 A) 
is only 4.7% smaller than the equilibrium lattice param- 
eter 6.013 A of RS GeTe (both values quoted are from 
calculations). Thus RS GeTe should act as a suitable 
substrate for growing RS CrTe films. 

In FigsHand[3]we show the DOSs for the RS and NA 
structure CrTe at their respective equilibrium volumes 
obtained via WIEN2k FP-LAPW code. Note that at 
equilibrium volume CrTe in NA structure is not half- 
metallic, but becomes so at expanded volumes. In FigJH 
the DOS of RS CrTe at the equilibrium RS-GeTe lattice 
parameter (6.013 A) is given, where it is seen to be half- 
metallic. The minority-spin gap is 1.388 eV, HM gap 
is 0.637 eV. Both the minority-spin gap and HM gap 
become larger with increasing lattice parameter and half- 
metallicity is robust with respect to lattice expansion. 

We have supplemented the FP-LAPW calculations 
with LMTO-ASA calculations using the tight-binding 
linear muffin tin orbital (TB-LMTO)^ method and 
its implementation for random alloys using the coher- 
ent potential approximation (CPA) i 51 ' 52 The exchange- 
correlation potential given by Vosko, Wilk and Nusair— 



TABLE I: Results obtained via the WIEK2A; code: equilibrium lattice constant a, Cr-Te bond- length Lcv-Te, cohesive energy 
E c and bulk modulus B for ZB, RS and NA CrTe in the FM state. E c is defined as negative of the total energy of the formula 
unit with respect to the sum of the energies of the free Cr and Te atoms. Ei, E2 and E m are as defined in the text. All results 
shown are obtained using GGA. In some cases the corresponding LDA values are shown in parentheses. 



Structure 


a (A) 


LcVTe(A) 


B(GPa) 


E c (eV) 


Ei(eV) 


E 2 (eV) 


E m (eV) 


ZB 


6.263 (6.088) 


2.712 


45.278 (57.962) 


6.13 


4.416 


0.284 


0.328 


RS 


5.727 (5.533) 


2.864 


59.744 (70.823) 


6.38 


1.605 


0.019 


0.076 


NA (c/a=1.542) 


4.118 


2.859 


46.487 


6.46 


1.453 


0.107 
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FIG. 2: (Color online) DOS for the RS structure at the equi- 
librium volume (a=5.7274 A) obtained via the FP-LAPW 
method. 
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FIG. 4: DOS in RS CrTe at the equilibrium lattice parameter 
of RS-GeTe (6.013 A). 
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FIG. 3: (Color online) DOS for the NiAs structure at the 
(a) equilibrium volume Vo (c/a=1.5424, a=4.118 A) and (b) 
10% expanded volume V=l.lVo (c/a= 1.5424), obtained via 
the FP-LAPW method. 



is used. We optimize the ASA (atomic sphere approx- 
imation) errors by including empty spheres in the unit 
cell. The DOS for the RS structure obtained via the 
LMTO-ASA method for three different lattice parame- 
ters, from just below to just above the critical value for 
the half-metallic phase is shown in FigfS] 

As pointed out above, several features of the electronic 
structure are dependent on the exchange-correlation ef- 
fects. Enhanced correlation generally increases the equi- 
librium lattice parameter and the magnetic moment for 
a given volume. The equilibrium lattice parameter for 
the ZB structure is 6.26 A according to the FP-LAPW 
calculations using GGA, 6.09 A according to FP-LAPW 
and LDA, and 6.03 A according LMTO-ASA method us- 
ing LDA. These LDA values value are in close agreement 
with the result 6.07 A from the Korringa-Kohn-Rostoker 
(KKR) Green function method and LDA M For RS CrTe 
equilibrium lattice parameter from FP-LAPW GGA is 
5.727 A, while replacing GGA by LDA results in the value 
5.53 A, a decrease of about 3%. Since GGA is not in- 
corporated at this stage in the LMTO-ASA code that we 
used, we checked the effect of enhancing the correlation 
effect in LMTO-ASA by using the LDA+U method^ 
The equilibrium lattice parameter for RS CrTe increases 



III. MAGNETIC MOMENTS 
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FIG. 5: LMTO-ASA DOS for the RS structure for three 
different lattice parameters, from just below to just above the 
critical value for the HM phase. 
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Results of the WIEN2/c calculations for magnetic mo- 
ments and some other useful quantities are given in Table 
HI1 The total magnetic moment contains contributions 
from the Cr atom, the Te atom and the interstitial region. 
The Cr atom carries the largest amount of magnetic mo- 
ment, the Te site has a small magnetic moment, about 
3-4% of that of the Cr atom and a certain amount of 
moment, 16-17% of the Cr-moment, is also found in the 
interstitial region. The magnetic moment in the intersti- 
tial region arises mainly from the d-states of the Cr atom. 
The moments of Cr and Te are antiparallel to each other, 
while the moments of the interstitial region and the Cr 
atom are parallel. The calculated total magnetic moment 
for both ZB and RS structures are 4.000 \ib per formula 
unit at their equilibrium lattice parameters. An integer 
value of the magnetic moment is a characteristic feature 
of HM ferromagnets. The ratio between the induced mo- 
ment on Te and the moment on Cr is 0.044, 0.031 and 
0.030 for ZB, RS and NA, respectively. The magnitudes 
of the moments on the two sublattices increase with in- 
creasing lattice parameters, due to decreased hybridiza- 
tion between Cr-d and Te-sp orbitals. Above a critical 
value of the lattice parameter, the moment per formula 
unit (f.u.) saturates at a value of 4.000 /is, as the HM 
state is achieved, while the local moments on the Cr- 
and Te-sublattices increase in magnitude, remaining op- 
posite to each other. The saturation value satisfies the 
so-called "rule of 8"iZ: M = (Z tot -8)fjt B , where Z tot is 
the total number of valence electrons in the unit cell. The 
number 8 accounts for the fact that in the HM state the 
bonding p bands are full, accommodating 6 electrons and 
so is the low-lying band formed of the s electrons from 
the sp atom, accommodating 2 electrons. The magnetic 
moment comes from the remaining electrons filling the 
d states, first the e g states and then the t2 g . The mo- 
ment versus volume variations for the three structures 
obtained in the WIEN2/c calculations are shown in Fig|7] 



FIG. 6: LMTO-ASA DOS for the RS structure using the 
LDA and LDA+U (U=1.02 eV) method for a lattice param- 
eter of 5.9423 A. 



from 5.245 A to 5.4% A, i.e. by about 3%, from LDA 
to LDA+U with U=1.02 eV (0.075 Ry) and J=0. Note 
that the width of the energy gap in the minority states for 
the RS structure also has a noticeable dependence on the 
treatment of exchange and correlation effects. This gap 
is widened with enhancement of correlation effects. For 
example, for a lattice parameter just above the critical 
value for the half- metallic state (Fig®, the LMTO-ASA 
calculation with LDA gives a gap width of approximately 
0.5 eV, while LDA+U calculation for U=0.075 Ry (1.02 
eV), the gap increases to ~0.8 eV, which compares well 
with the FP-LAPW GGA value of 0.76 eV. 



In LMTO-ASA method the interstitial region is de- 
scribed by empty spheres and the electronic properties 
of this region is primarily dictated by the tails of the or- 
bitals centered at the real atoms such as Cr and Te. As 
in the FP-LAPW calculation, we find that the magnetic 
moments of the empty spheres are parallel to those of 
the Cr atoms, while those of the Te atoms are antiparal- 
lel. For the RS structure, the moment on the Te atom is 
small (around 1-2% of that of the Cr atom) for low values 
of the lattice parameter, while the moment on the two 
empty spheres are each about 3.5-4% of that of the Cr 
atom. With increasing lattice parameter, the Te-moment 
increases with respect to the Cr-moment to about 7-8% 
at the highest lattice parameter studied. At these high 
lattice parameters the moments on the empty spheres, 
originating from the tails of the nearby Cr-3d orbitals, 
drop to about 2% of the Cr-moment. 

In the ZB structure, the two empty spheres are not 



TABLE II: Comparison of magnetic moments per formula unit (pb), minority-spin gap Grain (eV) and HM gap Ghm (eV) 
for FM CrTe in ZB, RS and NA structures, at their respective equilibrium volumes, obtained by using the FP-LAPW method 
(WIEN2fc). The numbers without and within the parentheses refer to values obtained by using GGA and LDA, respectively. 



Structure 


Cr 


Te 


Interstitial 


Total 


vimin 


Ghm 


ZB 

RS 

NA 


3.542 (3.501) 
3.563 (3.517) 
3.489 (3.398) 


-0.157 (-0.119) 
-0.112 (-0.079) 
-0.105 (-0.076) 


0.616 (0.619) 
0.548 (0.548) 
0.520 (0.503) 


4.000 (4.000) 
4.000 (3.986) 
3.904 (3.825) 


2.911 
0.761 


0.829 
0.028 
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FIG. 7: (Color online) Magnetic moment per formula unit 
(f.u.) for the NA, ZB and RS structures obtained via FP- 
LAPW (WIEN2fc) calculations. 



FIG. 8: (Color online) Magnetic moment per formula unit 
(f.u.) and of the Cr atoms in the ZB ans RS structures ob- 
tained in the LMTO-ASA calculations. 



equivalent and the induced moments on the empty 
spheres and the Te atoms in this case has been described 
in detail in our recent publication.— For the ZB struc- 
ture the induced moment on Te is negligible for small 
lattice parameters, increases steadily with increasing lat- 
tice parameter, achieving a maximum value of 4% of the 
Cr-moment at the highest value of the lattice parameter 
studied. Most often it is less than 1%. The magnitudes 
of the induced moments on the empty spheres decrease 
in the following order: sublattice ES-1 (the sublattice 
of empty spheres that is at the same distance with re- 
spect to the Cr-sublattice as the Te-sublattice ), sublat- 
tice ES-2 (sublattice of empty spheres further away from 
the Cr-sublattice). The induced moment on ES-2 is al- 
ways small, decreases with increasing lattice parameter 
from about 2% to about 0.4% of the Cr-moment. The 
induced moment on ES-1 is higher than the other two in- 
duced moments for smaller lattice parameters. For larger 
lattice parameters it is comparable to the Te-moment, 
but opposite in sign. Above a critical value of the lattice 
parameter, the moment per formula unit (f.u.) saturates 
at a value of 4.0 /is, as the HM state is achieved, while 
the local moments on the Cr- and Te-sublattices increase 
in magnitude, remaining opposite in sign. Sandratskii 
et al^ have discussed the problem associated with such 



'induced moments' in calculating the Curie temperature. 
LMTO-ASA results using LDA for the ZB and RS struc- 
tures are presented in FigfE] 

Magnetism is strongly dependent on correlation effects. 
Even in the metallic states, where magnetism is believed 
to be described reasonably well by density functional 
methods and LDA, different treatment of the exchange- 
correlation effects lead to different values of magnetic mo- 
ments. We thus compare the magnetic moments of our 
LMTO-ASA calculations for LDA and LDA+U^i The 
enhanced correlation effects in the LDA+U calculations 
(using U=0.075 Ry, J=0) increase the values magnetic 
moments and bring the results in general close agree- 
ment with those obtained via using GGA in the WIE2fc 
package. This is shown in FigO The agreement between 
GGA and LDA+U results is due to the fact that the GGA 
goes beyond the LDA by including electron correlations 
in some way. 



IV. EXCHANGE INTERACTION AND CURIE 
TEMPERATURE 



Exchange interactions are usually calculated by 
mappin g 55 ' 56 the system energy onto a classical Heisen- 
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FIG. 9: (Color online) Magnetic moment per formula unit 
(f.u.) versus lattice parameter for RS CrTe calculated by us- 
ing LMTO-ASA method and LDA and LDA+U (U=1.02 eV, 
J=0) schemes (upper panel) . Enhanced correlation effects in- 
corporated via the LDA+U method brings the results in close 
agreement with those obtained via FP-LAPW-GGA method 
(lower panel). 



berg model: 



H c ff 



/ j Jij e i 



(1) 



where i,j are site indices, e^ is the unit vector pointing 
along the direction of the local magnetic moment at site i, 
and Jij is the exchange interaction between the moments 
at sites i and j. In this work, we use the method of Liecht- 
enstein et aL,— ~— based on multiple-scattering formal- 
ism, which was later implemented for random magnetic 
systems by Turek et al. , using CPA and the TB-LMTO 
methodi^i For a discussion of the various methods used 
in the mapping procedure the reader may consult our 
previous publication^ The exchange integral in Eq.([T]) 
is given by 



Jij = - — lim Im j trL 



4tt 



>o+ 



A l (z)gUz)A,(z)gi] dz , (2) 



where z = E + it represents the complex energy vari- 
able, L — (/,m), and Aj(z) = Pj(z) — Ff(z), repre- 
senting the difference in the potential functions for the 
up and down spin electrons at site Y. In the present 
work gfj{z)(a =t)4-) represents the matrix elements of 
the Green's function of the medium for the up and down 
spin electrons. The negative sign in Eq.([T]) implies that 
positive and negative values of Jij are to be interpreted 
as representing ferromagnetic and antiferromagnetic in- 
teractions, respectively. 

In the mean- field approximation (MFA), the Curie 
temperature T c , based on the exchange interaction be- 
tween the strong moments only, i.e. the Cr atoms, can 



FIG. 10: (Color online) Exchange interaction between the Cr 
atoms in the ZB and RS structures as a function of the lattice 
parameter. The two lattice parameters for the ZB structure 
chosen are above and below the equilibrium lattice parameter 
5.727 A for the RS structure, obtained in the FP-LAPW 
(WIEN2k) calculation. 



be written as. 



k R T 



AM'A _ 2 y^ jCr,Cr 



(3) 



where the sum extends over all the neighboring shells 
centered around the Cr-site 'o'. An improved descrip- 
tion of finite-temperature magnetism is provided by the 
random phase approximation (RPA), 62 with T c given by 



[ki 



T^r 1 = ~ £[ jCr ' Gr (°)- J * ^]- 1 ■ (4) 



Here N denotes the order of the translational group ap- 
plied and J ' (q) is the lattice Fourier transform of 
the real-space exchange integrals J i r ' r . T c iPA can be 
shown to be smaller than T^ AFA ^- and has been found to 
be in close agreement with results given by Monte- Carlo 
simulations. 63 Unless otherwise stated, the T c estimates 
given in this paper are based on RPA, as given by Eq.(|l|). 
Sandratski: 55 et al. have discussed the case when, in 
addition to the interaction between the strong moments, 
there is one secondary, but much weaker, interaction be- 
tween the strong and one induced moment. In this case, 
the Curie temperature, calculated under the mean-field 
approximation (MFA) , is enhanced due to this secondary 
interaction, irrespective of the sign of the secondary in- 
teraction. In other words, the Curie temperature would 
be somewhat higher than that calculated by consider- 
ing only the interaction between the strong moments. 
The corresponding results under the random phase ap- 
proximation (RPA) have to be obtained by solving two 
equations simultaneously. One can assume that the RPA 
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results for the Curie temperature follow the trends repre- 
sented by the MFA results, being only somewhat smaller, 
as observed in the absence of induced moments. 

In Fig llOl we compare the exchange interactions be- 
tween the Cr atoms obtained by using Eq.© for the equi- 
librium lattice parameter of 5.727 A (FP-LAPW GGA 
result) of RS CrTe with those in ZB CrTe for two values 
of the lattice parameter, one above and the other below 
the value 5.727 A. For the RS structure the interactions 
are all ferromagnetic, while for the ZB structure they 
are predominantly antiferromagnetic in this range of the 
lattice parameter. In our previous publication we have 
shown that due to these antiferromagnetic interactions, 
the Curie temperature in ZB CrTe is drastically reduced 
for a certain range of the lattice parameter. The ground 
state for this structure may even have a more complex 
magnetic phase than simple ferromagnetic. Thus for lat- 
tice parameters around 5.74-5.75 A RS CrTe is not only 
guaranteed to be a half-metallic ferromagnet, but also 
should have a higher T c than ZB CrTe of similar volume 
per atom. 

In FigfTT]we compare the Tf- PA values for RS and ZB 
CrTe obtained by using LMTO-ASA and LDA. The ref- 
erence state used for the calculation is the FM state. 
For the ZB phase, the use of Eq.(4) neglecting all the 
induced moments, results in unphysical negative T c in 
the lattice parameter range ~5. 7-6.0 A because of dom- 
inant antiferromagnetic Cr-Cr exchange at the nearest 
and next nearest neighbor sites. The ZB phase has been 
discussed in detail in our previous publication^ where 
we showed that in spite of these nearest and next near- 
est neighbor AFM interactions, the FM state is lower in 
energy than the AFM[001] and AFM[111] states. The 
cumulative effect of the FM interactions involving the 
more distant pairs of Cr atoms, as well as that due to 
the induced moments, may stabilize the FM phase. Ei- 
ther the cumulative effects of distant pair interactions 
make the system FM or the magnetic state is other than 
AFM[001] and AFM[111], possibly complex. If the mag- 
netic state is indeed FM, then according to the prescrip- 
tion of Sandratskii^ et al. T c should be enhanced with 
respect to the value given by Eq.(4) and may be positive 
but small. This explains why Sreenivasan et alM and 
Bi et oi4S observe T c s only as high as 100 K in their 
thin film sample of ZB CrTe. For RS CrTe, on the other 
hand, all dominant interactions are FM, and as a re- 
sult the T c calculated by using Eq.(4) and based on only 
Cr-Cr interactions give positive values. The actual T c s 
are most probably higher than what is shown in Fig llll 
The effect of enhanced correlation for RS CrTe incorpo- 
rated via LDA+U method, using U=0.075 Ry, and J=0, 
is shown in Fig |12l Within the uncertainties of our re- 
sults there seems to be no appreciable difference between 
the predictions of LDA and LDA+U. 

It would transpire from FigfTU] that a J1-J2 model 
may suffice to describe the magnetism of the RS phase. 
However, for quantitative accuracy of T c the more dis- 
tant interactions, although individually small compared 



to the first two, cannot be neglected, i.e. their cumu- 
lative effect is non-negligible. Note that for the ZB 
phase a J1-J2 model would clearly be inadequate. For 
the RS phase, if one is prepared to neglect interactions 
beyond the next nearest neighbor pairs, then the near- 
est neighbor interaction Ji and the next nearest inter- 
action J2 can be obtained from the energies of the FM 
and AFM[001] and AFM[111] states^ We have calcu- 
lated Ji and J2 using this procedure for the FP-LAPW- 
GGA equilibrium lattice parameter of 5.727 A, and the 
results are: Ji=3.76 meV, J 2 =27.63 meV (FP-LAPW- 
GGA), Ji=3.29 meV, J 2 =18.32 meV (FP-LAPW-LDA), 
Ji=4.61 meV, J 2 =20.0 meV (LMTO-ASA-LDA). To 
compare these results with those given by the procedure 
based on the multiple scattering formalism and Eq.Q we 
need to consider that all the distant neighbor interactions 
between the Cr atoms (and the interactions involving the 
induced moments) get mapped effectively on to just two 
interactions in the J1-J2 model. This is tantamount to 
reproducing the band structure of a long-ranged Hamilto- 
nian via a TB-model with only nearest and next-nearest 
neighbor interactions. Success is possible, but not always 
guaranteed (unless further interactions are included). A 
meaningful way to compare the two approaches is to 
compute the dispersion J(q) = 5Z q </ r r exp(iq-R). 
In FigfTni we show these dispersion curves for RS CrTe 
along the various symmetry directions and all the dif- 
ferent cases considered. The values of the exchange in- 
teractions quoted above for the J1-J2 model would put 
the T c for the lattice parameter 5.727 A around ^1000 
K. One reason that the magnitudes of the J(q) values 
based on Eq.Q come out smaller than those given by 
the J1-J2 model is that in computing J(q) from Eq.([2|) 
we have neglected all interactions involving the induced 
moments. 

In our previous publication^ we have argued that esti- 
mates of T c from above the transition can be obtained by 
considering the reference state for the exchange calcula- 
tion in the method of Liechtenstein et alM-~— to be the 
paramagnetic state described by the disordered local mo- 
ments (DLM)££~— model. In this description magnetic 
moments appear only on the magnetic, i.e. Cr, atoms, 
but are oppositely directed for half of the randomly cho- 
sen atoms, giving a zero moment in total. The random- 
ness is treated using the CPA^i^ T c given by the DLM 
reference states are usually higher than their FM coun- 
terparts. The DLM estimates of T c are free from errors 
originating from the neglect of the exchange interactions 
involving induced moments in Eqs.([3]) and ^. For this 
reason we have computed the T c using the DLM model as 
well and compare these with the T c obtained for the FM 
reference states in Figfn] As expected, the DLM val- 
ues are usually higher for most of the lattice parameters. 
Since the ground state of RS CrTe is known to be FM 
for the entire range of the lattice parameters studied, we 
expect the FM results to be valid, although somewhat un- 
derestimated due to the neglect of the induced moments. 
Thus the theoretical estimates of T r should lie between 
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FIG. 11: (Color online) The RPA Curie temperature (Eq.©) 
as a function of the lattice parameter for the ZB and RS struc- 
tures of CrTe, calculated by considering the Cr-Cr interaction 
only and by using the FM reference state. The arrows indicate 
the equilibrium lattice parameter values according to various 
calculations. The results for T c are somewhat affected by 
the 'induced' moments on the Te atoms as well as the empty 
spheres used in the LMTO-ASA calculation. For ZB CrTe, 
this error, combined with strong antiferromagnetic spin fluc- 
tuation for the lattice parameters in the range ~5. 7-6.0 A, 
result in unphysical negative values of T c in this range (see 
text for details). Away from this range T c values are reliable, 
although somewhat underestimated. 
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FIG. 12: (Color online) Comparison of the Curie temperature 
in RS CrTe calculated for the FM reference states in RPA us- 
ing the LMTO-ASA scheme and LDA and LDA+U methods. 



FIG. 13: (Color online) Comparison of the lattice Fourier 
transforms of the Cr-Cr exchange interactions in RS CrTe for 
the FP-LAPW-GGA equilibrium lattice parameter of 5.727 
A and the various cases discussed in the text. 



those given by the FM and DLM results. Experimen- 
tal results on Cr-based chalcogenides and pnictides show 
that these estimates are on the high side. One should 
note, however, that second order phase transition in a 
bulk material is different from that in thin films grown on 
substrates and T c s may differ even though the two tran- 
sitions belong to the same universality class. A point to 
note is that for larger volumes per atom the DLM energy 
is lower than the NM energy, with a Cr-moment that 
is slightly less than that for the FM state of the same 
volume. Dominating J^-s for the DLM reference states 
at these volumes are AFM-like. For smaller volumes the 
Cr-moment disappears, indicating strong hybridization 
between the orbitals on Cr- and Te-sublattices. 



V. SUMMARY OF RESULTS 

We have provided a comparison of the magnetic prop- 
erties of CrTe in three different structures, NA, RS and 
ZB, at their respective equilibrium volumes. A detailed 
comparison of the half-metallic properties of the ZB and 
RS phases is provided, as thin films grown on most sub- 
strates are likely to have either of these two structures. 
The RS phase is more stable than the ZB phase, in that 
the lowest energy state in the RS phase is about 0.25 
eV/atom lower than that for the ZB phase. For some 
range of the lattice parameter, namely between 5.6-6.1 
A , in the ZB phase there are significant antiferromag- 
netic interactions between the Cr atoms. The RS phase 
appears to be free from such strong antiferromagnetic in- 
teractions. Ferromagnetism appears to be more robust 
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FIG. 14: (Color online) Variation T c as a function of lattice 
parameter for FM and DLM reference states in RS CrTe 



in the RS phase than in the ZB phase. Increased antifer- 
romagnetic interactions in the ZB phase result from the 
smaller distance between the Cr and the Te atoms and 
consequently increased hybridization between the Cr-3d 
and Te-5p orbitals in the ZB phase than in the RS phase 
for the same lattice parameter. Increased stability with 
respect to the ZB phase, stronger ferromagnetism and 
the promise of a higher T c should make the growth and 
measurements of the physical properties of RS CrTe an 
attractive project for experimentalists in the field. 
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